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■ Abstract. The time delays between the components of a lensed quasar are basic tools to 

r/3 . 

analyze the expansion of the Universe and the structure of the main lens galaxy halo. In this 
paper, we focus on the variability and time delay of the double system SBS 0909+532A,B 
as well as the time behaviour of the field stars. We use VR optical observations of SBS 
0909+532A,B and the field stars in 2003. The frames were taken at Calar Alto, Maidanak 
and Wise observatories, and the VR light curves of the field stars and quasar components are 
derived from aperture and point-spread function fitting methods. We measure the /?-band 
time delay of the system from the x 1 and dispersion techniques and 1000 synthetic light 
curves based on the observed records. One nearby field star (SBS 0909+532c) is found to 
be variable, and the other two nearby field stars are non-variable sources. With respect to 
the quasar components, the R-hanA records seem more reliable and are more densely pop- 
ulated than the l/-band ones. The observed ^-band fluctuations permit a pre-conditioned 
measurement of the time delay. From the x 2 minimization, if we assume that the quasar 
emission is observed first in B and afterwards in A (in agreement with basic observations 
of the system and the corresponding predictions), we obtain At ba = — 45 * l n days (95% 
confidence interval). The dispersion technique leads to a similar delay range. A by-product 
of the analysis is the determination of a totally corrected flux ratio in the R band (corrected 
by the time delay and the contamination due to the galaxy light). Our 95% measurement 
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Amba = m B(t + At sa ) - m A (t) = 0.575 ± 0.014 mag is in excellent agreement with previous 
results from contaminated fluxes at the same time of observation. 

Key words. Gravitational lensing - Quasars: general - Quasars: SBS 0909+532 - Stars: 
variables: general 

1. Introduction 

The system SBS 0909+532 was discovered by Stepanyan et al. (1991). Some years later, 
a collaboration between the Hamburger Sternwarte and the Harvard-Smithsonian Center for 
Astrophysics resolved the system into a pair of quasars (A and B) with a direct 7?-band flux 
ratio (at the same time of observation) Am = nig — iha = 0.58 mag and a separation of about l'.'l 
(Kochanek et al. 1997). The direct /?-band flux ratio was not consistent with the direct flux ratios 
at other wavelengths: Am = 0.31 mag in the / band and Am = 1.29 mag in the B band. From 
observations with the 4.2 m William Herschel Telescope, a Spanish collaboration got spectra for 
each component of the system. The data showed that the system consists of two quasars with the 
same redshift (z s = 1 .377) and identical spectral distribution, supporting the gravitational lens in- 
terpretation of SBS 0909+532 (Oscoz et al. 1997). Oscoz et al. (1997) detected a Mgn doublet in 
absorption at the same redshift (z a bs = 0.83) in both components, and they suggested that the ab- 
sorption features were associated with the photometrically unidentified lensing galaxy. Through 
a singular isothermal sphere (SIS) lens model, the authors also inferred the first constraint on 
the time delay between the components: |AtbaI < 140 days, where Atba is the delay of B with 
respect to A and the Hubble constant is assumed to be Hq = 70 km s _I Mpc -1 . 

In recent years, Lubin et al. (2000) indicated the possible nature of the main deflector (early- 
type galaxy) and confirmed its redshift (zd = 0.830). Lehar et al. (2000) reported on a program 
including Hubble Space Telescope (HST) observations of SBS 0909+532. They discovered the 
main lens galaxy between the components, which has a large effective radius, with a correspond- 
ingly low surface brightness. This lens galaxy is closer to the brightest component (A), which 
is not in contradiction with SIS-like lens models when the farther and fainter component (B) is 
stronger affected by dust extinction (see below). The colors of the lens are consistent with those 
of an early-type galaxy at redshift 0.83. Assuming a singular isothermal ellipsoid (SIE) model, 
Lehar et al. predicted a time delay At B a in the range [- 10, - 87] days (Hq = 70 km s -1 Mpc -1 ). 
At a given emission time, the sign "-" means that the corresponding signal is observed first in 
B and later in A. The COSMOGRAIL collaboration provided the distribution of predicted time 
delays of the system (Saha et al. 2005). In their histogram (Fig. 10 of Saha et al.), there are two 
features: the main feature is an asymmetric peak around - 80 days and the secondary one is 
another asymmetric peak around - 45 days. Therefore, if the COSMOGRAIL predictions are 
right, the time delay is very probably of 2-3 months (component B leading component A), but 
we cannot rule out a delay of about one and a half months. On the other hand, the flux ratio 
Send offprint requests to: A. Ullan 
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anomaly pointed by Kochanek et al. (1997) was confirmed and accurately studied by Motta et 
al. (2002) and Mediavilla et al. (2005), who reported the existence of differential extinction in 
the main lens galaxy. Chartas (2000) and Page et al. (2004) also studied the system in the X-ray 
domain. 

Time delays are basic tools to discuss the present expansion rate of the Universe and the struc- 
ture of the main lens galaxy haloes (e.g., Refsdal 1964; Kochanek, Schneider & Wambsganss 
2004), so that variability studies are crucial. While some time delays have been measured 
from radio light curves (PKS 1830-211: Lovell et al. 1998; Q0957+561: Haarsma et al. 1999; 
B0218+357: Biggs et al. 1999; B1600+434: Koopmans et al. 2000; B1422+231: Patnaik & 
Narasimha 2001; B 1608+656: Fassnacht et al. 2002) or X-ray variability (e.g., Q2237+0305: 
Dai et al. 2003), an important set of delays are based on optical monitoring of gravitation- 
ally lensed quasars. Optical frames taken at Apache Point Observatory, Fred Lawrence Whipple 
Observatory and Teide Observatory were used to estimate a 14-month delay for the double sys- 
tem Q0957+561 (e.g., Pelt et al. 1996; Kundic et al. 1997; Serra-Ricart et al. 1999; Ovaldsen et 
al. 2003). Although the time delay of this first multiple quasar has been confirmed through inde- 
pendent observations, the measurement is only 5% accurate, or equivalently, there is an uncer- 
tainty of about 20 days (Goicoechea 2002). The Tel- Aviv University (TAU) group have recently 
determined the time delay between the two components of HE 1 104-1805 (Ofek & Maoz 2003). 
The TAU delay of HE 1 104-1805 disagrees with the earlier estimation by Gil-Merino, Wisotzki 
& Wambsganss (2002), but it is in excellent agreement with the determination by Wyrzykowski 
et al. (2003). Schechter et al. (1997) measured two delays for the quadruply imaged quasar PG 
1 1 15+080. The Belgian-Nordic collaboration carried out a very intense activity during the past 
five years. They participated in several monitoring projects and measured several time delays at 
optical wavelengths: B1600+434 (Burud et al. 2000), HE 2149-2745 (Burud et al. 2002a), RXJ 
0911.4+0551 (Hjorth et al. 2002), SBS 1520+530 (Burud et al. 2002b) and FBQ 0951+2635 
(Jakobsson et al. 2005). The formal accuracies of these 5 estimations range from 5 to 25% (the 
lcr error bars vary from 4 to 24 days). Kochanek et al. (2005) also measured the time delays 
between the components of the quadruple quasar HE 0435-1223. 

The aim of this paper is to present VR observations of SBS 0909+532 in 2003 conducted 
by the University of Cantabria (UC, Spain), the Institute of Astronomy of Kharkov National 
University (lAKhNU, Ukraine), the Sternberg Astronomical Institute (SAI, Russia) and the Ulug 
Beg Astronomical Institute of Uzbek Academy of Science (UBAI, Uzbekistan). We also present 
TAU observations of the field stars in 2003, which have been kindly made available to us. This 
new optical monitoring campaign was carried out at the Calar Alto Observatory (Spain), the 
Maidanak Observatory (Uzbekistan) and the Wise Observatory (Israel), and the frames were 
taken with the 1.5 m Spanish telescope, the 1.5 m AZT-22 telescope at Mt. Maidanak and the 
Wise Observatory 1 m telescope (Section 2). In Section 3, we describe the methodology to obtain 
the fluxes of the quasar components and the field stars. The VR light curves are also shown in 
Section 3. Section 4 is devoted to the time delay estimation from the light curves of A and B 
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Fig. 1. Calar Alto image of SBS 0909+532, showing a FOV of ~ 7' x 7'. The field contains the 
gravitationally lensed quasar ("GL") and six bright and non-saturated stars ("a-c", "sl-s2" and 
"x"). The nearby ("a-c") and relatively far ("sl-s2") field stars were introduced by Kochanek et 
al. (1997) and Nakos et al. (2003), respectively. A sixth field star ("x") is also included in the 
FOV. 

(quasar components) in the R band. Finally, in Section 5 we summarize our conclusions and 
discuss the feasibility of an accurate determination of the cosmic expansion rate and the surface 
density in the main lensing galaxy. 

2. Observations 

We have three different sets of frames for SBS 0909+532. The first set of optical frames cover 
the period between 2003 March 4 and June 2, and they are part of a UC project to test the 
feasibility of quasar monitoring programs through 1 -2 m telescopes in Spain (Ullan 2005). These 
observations were made with the 1 .52 m Spanish telescope at Calar Alto Observatory (EOCA), 
Almeria, Spain (see Ziad et al. 2005 for a site-testing on Calar Alto). The EOCA is equipped 
with a Tektronics 1024x1024 CCD detector, which has pixels with a physical size of 24 //m, 
giving a 0.4 arcsec pixeF 1 angular scale. The gain is 6.55 e~/ADU and the readout noise is 6.384 
e~. During this first monitoring, exposures in the V and R Johnson-Cousins filters were taken 
every night when clear, what makes a total of 20 observing nights. Bad weather in 2003 March 
and April prevented us from achieving a very dense sampling. For each monitoring night we have 
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Fig. 2. Maidanak zoomed-in image of SBS 0909+532 (about 8"on a side). We choose one of the 
frames taken with the best seeing conditions, and then we expand the portion occupied by the 
two quasar components A (the brightest object) and B (the faintest object), i.e., the GL in Fig. 1. 
To avoid the pixellation effect when a subframe is expanded, the brightness distribution is also 
smoothed. 

three consecutive frames on each filter, i.e., three 300 s exposures in the V passband and three 
180 s exposures in the R passband. Those were the maximum exposure times to avoid saturation 
of selected stars in the field. In Figure 1 we show a typical frame. In this typical exposure, half 
a dozen bright and non-saturated stars were fitted within the field of view (FOV). Following the 
notation of Kochanek et al. (1997), the FOV included the gravitationally lensed quasar ("GL") 
and nearby field stars "a" (South), "b" (North) and "c" (West). The FOV also included two stars 
that were introduced by Nakos et al. (2003) and were labelled as "si" and "s2". These two stars 
are placed relatively far from the gravitational lens system, and they appear close to the North- 
West edge of the frame (see Fig. 1). A sixth star ("x") appears close to the South-West edge of 
the typical frame. 

The second set of observations include frames in February 2003 as well as during April-May 
and October-November 2003. The total number of nights is 18. In this second program the im- 
ages were taken with the 1 .5 m AZT-22 telescope at Maidanak Observatory (Uzbekistan), with 
near diffraction-limited optics and careful thermo-stabilization, which allow for high-angular- 
resolution imaging. The AZT-22 telescope has a LN-cooled (liquid nitrogen cooled) CCD- 
camera, SITe- 005 CCD, manufactured in Copenhagen (Denmark). For this camera, the imaging 
area is split into 2000x800 pixels, where the pixel size is 15 fim and the intrinsic angular scale 
is 0.26 arcsec pixel -1 . The frames were taken in the R Bessel filter, which corresponds approx- 
imately to the R Johnson-Cousins passband. The poor tracking system of this telescope allows 
only exposures up to 3 minutes. To obtain sufficiently high photometric accuracy, we took sev- 
eral frames each observation night. With respect to the rectangular FOV of the telescope, the 
North/South coverage was 2.5 times smaller than the East/West one, so the "si", "s2" and "x" 
stars were not included within the FOV. Figure 2 shows a zoomed-in image made from one of 
the best frames in terms of seeing. There are two close quasar components, but the very faint 
galaxy is not apparent. The observations at Mt. Maidanak are part of IAKhNU, SAI and UBAI 
projects to follow up the variability of gravitationally lensed quasars. 



6 



A. Ullan et al.: Time delay of SBS 0909+532 



Observatory (Telescope) Frames/night (Filter) Observation Periods 

Calar Alto (1.5 m) 3 x 300 s ( V) + 3 x 180 s (R) March-June 

Maidanak (1.5 m) (3-11) x 180 s (R) February, April-May, October-November 

Wise (1.0 m) 1 x 420 s (R) 22 unevenly distributed nights 

Table 1. Observations of SBS 0909+532 in 2003. 



For the past six years the TAU group have been monitoring several gravitationally lensed 
quasars with the Wise Observatory 1 m telescope. The targets are mainly monitored in the 
Johnson-Cousins /?-band, and the frames are obtained with a cryogenically cooled Tektronix 
1024x 1024-pixel back-illuminated CCD. The angular scale is 0.7 arcsec per pixel. This pixel 
scale and the median seeing (FWHM) of about 2" do not allow resolving most of the lensed 
objects, e.g., SBS 0909+532. However, the frames of SBS 0909+532 in 2003 are characterized 
by wide FOVs, which incorporate the "a-c", "sl-s2" and "x" stars. This fact permits to do differ- 
ential photometry between several pairs of field stars, and thus, to test the reliability of the Calar 
Alto and Maidanak records. 

The pre-processing of the images included the usual bias subtraction, flat fielding using 
sky flats, sky subtraction and cosmic ray removal by using the Image Reduction and Analysis 
Facility (IRAF) and Munich Image Data Analysis System (MIDAS) environments. Some de- 
tails about the whole observational campaign are included in Table 1 (observatories, telescopes, 
frames/night, filters and observation periods). 

3. Photometry and VR light curves 

Due to the small angular separation between the two lensed components, about 1'.' 1 (Kochanek 
et al. 1997), the photometry of SBS 0909+532 is a difficult task. This task is also complicated 
by the presence of the main lensing galaxy between the components, which could make the 
computation of individual fluxes even harder. In general, aperture photometry does not work, 
so we must look for better approaches. An initial issue is to decide about the inclusion or not 
inclusion of a photometric model for the lensing galaxy. In principle, when computing the fluxes 
of SBS 0909+532 we may use a galaxy model derived from the HST images of the system. 
The galaxy model could also be inferred from the best images in terms of seeing. Once the 
relevant information on the galaxy is known, we would apply a PSF fitting method to all optical 
images, setting the galaxy properties to those derived from the HST or the best-quality images, 
and allowing the remaining parameters to vary (e.g., McLeod et al. 1998; Ullan et al. 2003). 
Magain, Courbin & Sohy (1998) also presented an alternative task (deconvolution) that combines 
all the frames obtained at different epochs to determine the numerical light distribution of the 
lensing galaxy as well as the positions of the point-like sources (quasar components), since 
these parameters do not vary with time. The flux of the point-like sources are allowed to vary 
from image to image, which produces the light curves. However, these and other procedures have 
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a reasonable limitation: they only work well when the galaxy light has a significant contribution 
to the crowded regions in the individual frames. For a very faint galaxy in a standard (i.e., not 
superb) frame, there is confusion between galaxy signal and noise, so the use of a given galaxy 
model could lead to biased fluxes of the components. The biases will depend on the quality of the 
image (seeing, signal-to-noise ratio, etc), which must produce artificial variability superposed 
to the real one. On the other hand, the use of a direct PSF fitting method (neglecting the galaxy 
brightness) leads to contaminated fluxes of the components. But if the galaxy is very faint, the 
contaminations will be small. Moreover, the variation of the quasar fluxes, seeing conditions, 
etc, will cause fluctuations in the contaminations, which are expected to be below the typical 
contamination levels. For standard frames of a quasar lensed by a very faint extended object, it 
is really difficult to choose between both approaches (with and without galaxy). 

Most of the Calar Alto and Maidanak individual frames of SBS 0909+532 do not show evi- 
dences for a galaxy brightness profile. This fact is due to the faintness of the galaxy, as we corrob- 
orate here below. If we consider an hypothetical astronomer that neglects the galaxy brightness 
and does direct PSF fitting (without taking into account the galaxy when doing the computation 
of the fluxes), it is possible to attain a rough estimation of the maximum contamination from the 
galaxy to the closest component A (at 074 from the centre of the deflector). We take into account 
the paper about 10 lens systems by Lehar et al. (2000), where, in Table 3, we can find the best 
available photometric and astrometric (HST) data of SBS 0909+532. The authors were able to 
trace the galaxy light in the H passband, by measuring its position and brightness. If we use the 
colors in the same table, we conclude that m ga i ~ 19 mag and niA ~ 16 mag in the / band (near- 
IR), and m ga i > 20.4 mag and tha ~ 16.7 mag in the V optical band. Therefore, as the R filter is 
placed just between the / filter and the V one, we may assume that m ga i - niA ~ 3.5 mag in the R 
band. The difference of 3.5 mag is consistent with a ratio of fluxes F ga i/FA of about 1/25. Thus, 
in the case of QSO 0957+561 we found a /?-band ratio of fluxes F ga i/FA of about 1/2.5 (Ullan et 
al. 2003), and now we have F ga i/FA ~ 1/25, what explains our unsuccessful efforts when mea- 
suring the flux of the lens galaxy in standard frames. As a result of that, in an extreme case (when 
direct PSF fitting leads to a magnitude niA+ ga i instead of niA, i.e., all the galaxy light is included 
in the profile of the A component) we find a relationship: niA = niA+ ga i + F ga i/FA, where the true 
flux (in magnitudes) niA differs from the contaminated flux through direct PSF fitting (mA+ ga i) in 
a quantity F ga [/FA- This maximum contamination of A would be only of 40 mmag, and the real 
contamination of both components will be less than our upper limit. The artificial fluctuations 
(caused by variable contamination) will be even smaller than the typical contamination levels, so 
we expect they will not play an important role in analyses of quasar variability (e.g., time delay 
estimates). 

In order to derive the light curves of the components A and B, we decide to use a direct 
PSF fitting method and do not consider the galaxy brightness in the fits. The key idea of this 
procedure is to obtain the different fluxes we are interested in by using a PSF that comes from 
a bright star in the field common to all frames. The point-like objects (quasar components and 



8 



A. Ullan et al.: Time delay of SBS 0909+532 




1.5 2 2.5 

FWHM (arcsec) 



to - 



B 2 h 



C/2 



O 



V filter 



A A A 
A 



I A 



• A. 



4 I A i 



A 
A 



1 1.5 2 2.5 3 

FWHM (arcsec) 

Fig. 3. (S /N) max versus FWHM for the Calar Alto frames in the R filter (top panel) and the V 
filter (bottom panel). The circles and triangles represent the frames with good and bad post-fit 
residues, respectively (see main text). 



stars) are modelled by means of the empirical PSF. Hence, we do not use a theoretical PSF (i.e., 
Gaussian distribution, Lorentzian distribution, etc), but the two-dimensional profile of a star in 
this field (a PSF star). Apart from a PSF star, we also need a reference star to do differential 
photometry and to obtain relative fluxes Ma - m re f and Mb - m re f- The good behaviour of the 
reference star is usually checked by using a control star, so the fluxes m con - m re f are expected 
to agree with a constant level. Nevertheless, since the /?-band flux ratio is discussed in Section 
5, we also want to obtain a rough estimation of the contaminations from this direct technique. 
With this aim, a deconvolution technique (Koptelovaet al. 2005) is also applied to a set of frames 
with good seeing and signal. The selected frames are fitted to a model including the galaxy, and 
thus, we are able to obtain a few clean fluxes of components A and B and compare them with the 
corresponding contaminated fluxes (through a direct PSF fitting). The averaged contaminations 
are used in Section 5. 
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Fig. 4. Calar Alto light curves of the double quasar in the R filter. We use the "b" star as the 
reference object, because it is confirmed as a non-variable star in Fig. 7. The circles are the 
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Fig. 5. Calar Alto light curves of the double quasar in the V filter. We use the non-variable star 
"b" as reference object (see Fig. 7). The circles are the fluxes yA and the squares are the fluxes 
ya— 0.65 mag. 



3.1. Calar Alto frames and light curves 

We adopt a model of the system including two point-like sources and a constant background. 
The model is fitted to each image by adjusting its 7 free parameters (two-dimensional positions 
of A and B, instrumental fluxes of both components and background) to minimize the sum of the 
square residuals, as described in McLeod et al. (1998) and Lehar et al. (2000). We use windows 
of 64x64 pixels. Each empirical PSF is a subframe of 64x64 pixels around the PSF star (the "a" 
star in Fig. 1), while the lens system is analyzed from a subframe of the same size, but centered 
on the double quasar. The instrumental fluxes of the "b", "c", "si" and "s2" stars are also inferred 
from 64x64 pixels windows centered on them. We initially focus on the nearby field stars, and 
take the "b-c" stars as the control-reference objects. The "a" object is the brightest star in the 
"a-c" triangle, and "b" and "c" were spectroscopically identified by Kochanek et al. (1997) and 
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Fig. 6. Calar Alto light curves of the stars in the R filter. The Calar Alto fluxes are checked 
from six Wise frames. Filled and open symbols are associated with PSF fitting and aperture, 
respectively. First, the non-variable star "b" (see Fig. 7) is taken as the reference object. The top 
box incorporates the behaviour of y„ +2.15 mag (y a — m a - nib): Calar Alto (open squares) and 
Wise (open triangles). The second box (under the top box) contains the fluxes y c (y c = m c - nib): 
Calar Alto (filled squares) and Wise (open circles). The third and fourth boxes include the y s \ + 
2.38 mag (y s \ = m s \ - nib) an d y$i + 0.29 mag (y S 2 = m s i - nib) records, respectively. In the third 
box: Calar Alto (filled and open triangles) and Wise (open astroids), whereas in the fourth box: 
Calar Alto (filled and open circles) and Wise (open rhombuses). Second, we plot the m S 2 - m s \- 
0.20 mag records in the bottom box: Calar Alto (filled and open star symbols) and Wise (open 
crosses). 



Zickgraf et al. (2003): "b" is a FG star, whose spectrum includes the G-band and Ca n H-K lines, 
and "c" is a M3 star. The R - I and B - R colors of the brightest component (A) and the "a" star 
are similar, the colors of the faintest component (B) are close to the colors of the "b" star, and 
the "c" star has colors different to those of the components and the "a-b" stars (see Table 1 of 
Kochanek et al. 1997). On the other hand, after checking the PSFs of the three nearby field stars 
("a-c"), we do not find significant differences between them. This suggests that the global shape 
of the PSFs around the lens system does not depend on the position and color of the point-like 
objects, so the PSF of the "a" star seems to be a reliable tracer of the PSF associated with any 
point-like object in the region of interest. 

As a first attempt for obtaining light curves we use the "b" and "c" stars as the control and 
reference objects, respectively. Unfortunately, we find clear evidences in favour of variability of 
the "c" star, since the three curves Ma - m c , Mb - m c and nib - m c have a similar global behaviour. 
This fact forces us to rule out the "c" star as a reference-control object and, thus, to take the 
"a" and "b" nearby field stars as the control and reference point-like sources, respectively. In the 
next subsection, we analyze the Maidanak-Wise fluxes m a - n% and show that both stars ("a" and 
"b") are non-variable objects. This result permits to assure the good behaviour of "b". The Calar 
Alto fluxes m a - 11% are not included in the analysis, since most Calar Alto data disagree with the 
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Maidanak-Wise common level of flux. We found an anomaly in the behaviour of the Calar Alto 
relative fluxes for widely separate stars (see below), so only the relative fluxes for neighbouring 
point-like objects are reliable photometric measurements. Fortunately, the comparison between 
the quasar components and the "b" nearby reference star seems to be a feasible approach. 

After applying the photometric method to the three individual frames for each filter and night 
(see Table 1), we obtain three different measurements of y A = m A - nib and y B = m B - nib in the 
V and R passbands for each night. To test the reliability of the instrumental fluxes of A and B, 
we analyse the residues in each residual frame. A residual frame is an image after subtracting the 
fitted background and point-like objects (PSF fitting method). More properly, we focus on the 
residual subframe occupied by the system, and then we estimate the residue-to-signal ratio (R/S ) 
in each pixel of interest. A R/S value less than 10% is acceptable, so a subframe with at least 90% 
of pixels having acceptable residues is considered to be related to reliable photometric solutions. 
Thus, we classify the individual fits in two categories: fits leading to < 90% of pixels having 
acceptable residues (bad fits, unreliable results) and good fits that are associated with reliable 
results (> 90% of pixels having acceptable residues). As a complementary test, we study the 
relation between the quality of the fits (in terms of post-fit residues) and two relevant parameters 
(image quality). The signal-to-noise at the brightest pixel of the lens system, (S /N) max , and the 
seeing, FWHM (in "), are the two parameters to compare with the fit quality. Some kind of 
correlation between good fits and good images is expected. In Figure 3 we draw the (S /N) max - 
FWHM plots for frames in the R filter (top panel) and the V filter (bottom panel). Circles and 
triangles represent good and bad fits, respectively. The plots in Fig. 3 indicate that the good fits 
correspond to images with high or moderate (S IN) max (> 30). Moreover, at moderate (S /N) max 
(~ 30-50), most of the good fits seem to be associated with a relatively good seeing (< 2 "). 
To obtain a robust photometry, we finally discard the frames corresponding to the triangles in 
Fig. 3. For each filter and night, if there are two or three good frames (good fits), then we get 
mean values of y A and y B , and compute standard deviation of means as errors. We only consider 
relative fluxes with uncertainties < 40 mmag. 

Now we plot y A (circles) and y B - 0.45 mag (squares) in Figure 4 (/?-band fluxes). If we 
concentrate our attention in the period with the best sampling (after day 2755), the A light curve 
shows a moderate decline and the B record shows a moderate rise. Indeed it seems that the 
"b" star is a good reference object (constant flux), since there is no zero-lag global correlation 
between y A and y B . In Figure 5 we show the light curves y A and y B - 0.65 mag in the V passband. 
In this case we have a total of 1 1 points for the A component (circles) and 10 points for the B 
component (squares). The V-band and /?-band light curves of the A component are consistent 
with each other. A final moderate decline appears in both curves. The situation is more confused 
for the B component. The /?-band final rise is not clearly reproduced in the V band, and the V- 
band final measurements could have underestimated formal errors. We note the relative faintness 
of B in the V band (Am ~ 0.8 mag), and thus, the possibility of systematic uncertainties when 
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the PSF fitting method is applied at some epochs. The data in both optical filters are available at 



After presenting the records of the double quasar, we concentrate on the Calar Alto light 
curves of the field stars that were previously introduced by Kochanek et al. (1997) and Nakos 
et al. (2003), i.e., y a — m a — nib, yc — tn c — nib, Js\ — m s\ — nib and y s i — m s i — nib- There 
are no previous studies on the variability of the nearby field stars "a-c". On the other hand, the 
farther field stars ("sl-s2") were verified to be non-variable by using 76 Wise frames taken from 
1999 December 24 to 2002 March 3 (Nakos et al. 2003). As Nakos et al. (2003) found that 
"si" and "s2" seem to be useful reference stars, we check the behaviour of "sl-s2" in 2003. 
The PSF of the stars in the surroundings of the double quasar could slightly differ from the PSF 
of the "sl-s2" stars in a relatively far region. Therefore, we must be careful when obtaining 
the instrumental fluxes of the farther stars. To detect possible anomalies caused by a mismatch 
between the brightness profile of the "a" star and the PSF of "sl-s2", the light curves y s i andy.,2 
are derived from both PSF fitting and aperture methods. The records y a , y c , y s \ and y.,2 in the R 
filter are depicted in Figure 6. To guide the eyes, we use some offsets and dashed horizontal lines 
and put all the relative records of each pair within a box. Filled and open symbols are associated 
with PSF fitting and aperture, respectively. The top box includes the y a + 2.15 mag fluxes (open 
squares). The second, third and fourth boxes (under the top one) correspond to the y c (filled 
squares), y s \ + 2.38 mag (filled and open triangles) and y S 2 + 0.29 mag (filled and open circles) 
records, respectively. As most of the stars are brighter than the quasar components (A and B) 
and they are far from other objects, the typical formal errors in the stellar fluxes are clearly less 
than the typical uncertainties in the fluxes of the components (these are usually fainter and are 
placed in a crowded region). The stellar error bars in Fig. 6 are often smaller than the sizes of the 
associated symbols. 

When doing aperture photometry on six 7?-band Wise frames covering the first semester 
of 2003, we obtain a y a + 2.15 mag light curve (open triangles in the top box of Fig. 6) that 
disagrees with the Calar Alto trend in the overlap period (between days 2710 and 2760). In the 
next subsection, we show that the Wise and Maidanak brightnesses are constant and consistent 
with each other, so the Calar Alto values of y a are not true fluxes, but anomalous results. On 
the contrary, the Wise light curve y c (open circles in the second box of Fig. 6) agrees with the 
Calar Alto curve in the overlap period. From the Wise frames we confirm the flux level during 
the high-state of "c". Unfortunately, the small-amplitude variability of "c" (rms fluctuation of ~ 
8 mmag) cannot be confirmed from the Wise data. The rms fluctuation of the Wise fluxes (~ 9 
mmag) is very similar to the Calar Alto variation, but the formal errors are relatively large (~ 10 
mmag). Moreover, there are no Wise frames in 2003 May (around the day 2780) and, thus, we 
cannot check (via Wise data) the reliability of the Calar Alto dip in y c (80-100 mmag). However, 
the flux of the "c" star at day 2793 in the V band confirms the existence of a transition from the 
low-state to the high-state, which is finished at days 2800-2810 (see the last open circle in the 
second box of Fig. 6). For the "sl-s2" stars, which are as far from star "b" as star "a" is, we 
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Fig. 7. Maidanak (open star symbols) and Wise (open triangles) light curves y a in the R filter 
(winter-spring of 2003). The solid line represents the global mean value and the two dashed 
lines describe the rms fluctuation of the measurements. This rms variation agrees with the typical 
uncertainty, so "a" and "b" seem to be non-variable sources. 



again find a disagreement between the Calar Alto trends and the Wise records (open astroids and 
rhombuses in the third and fourth boxes of Fig. 6). Although aperture curves are closer to the 
Wise behaviours, we cannot fairly reproduce the Wise data. Some probes with the "x" star (using 
y.x = tn x - Mb) also indicate that the Calar Alto and Wise behaviours disagree. It seems that the 
differential photometry between widely separate stars may lead to meaningless results, and only 
the relative fluxes for neighbouring objects are reliable. To test this conclusion, apart from the 
successful results through the neighbouring stars "b" and "c", we also analyze the differential 
photometry between the pair "sl-s2" (see Fig. 1). The curves m S 2 - m s \— 0.20 mag are depicted 
in the bottom box of Fig. 6: Calar Alto (filled and open star symbols) and Wise (open crosses). 
In the overlap period (from day 2710 to day 2760), there is a reasonable agreement between 
the results from both observatories, and the Calar Alto measurements seem to be quite reliable. 
From the Calar Alto frames, both photometric techniques are consistent with each other, but a 
constant flux cannot explain the observations. When we fit the data sets to a constant, our best 
solutions are characterized by x 1 ~ 162 (PSF fitting) and^ 2 ~ 6 (aperture). It is a curious fact 
that aperture photometry on only one frame per night leads to relative fluxes in rough agreement 
with a constant level. However, more refined measurements (aperture or PSF fitting on several 
frames per night) reveal the variability of one ("si" or "s2") or both stars. 

3.2. Maidanak frames and global R-band light curves of SBS 0909+532 

In the case of the R-band Maidanak observations, in order to derive the relative fluxes of the 
components of SBS 0909+532, we also use a direct PSF fitting. For a given frame, after to obtain 
a first estimate of the free parameters (initial solution), the fit is refined through an iterative 
procedure, which works as the CLEAN algorithm (0stensen 1994). The iterative task is done 
with each individual image, and the solutions converge after a few cycles. For each night, we take 
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Fig. 8. Global fl-band fluxes of SBS 0909+532 in 2003. The open circles (Maidanak) and filled 
circles (Calar Alto) are the fluxes y&, whereas the open squares (Maidanak) and filled squares 
(Calar Alto) are the relative fluxes yg- 0.45 mag. The top panel contains the results in the winter- 
spring of 2003 and the bottom panel includes the results in the autumn of 2003. 



all the available images and obtain the mean values of y^ and ys- From the standard deviation of 
the means, we also derive the errors in yA and ys- In agreement with the criteria in subsection 3.1, 
only fluxes with errors less than or equal to 40 mmag are considered. Apart from the analysis 
of the lens system, using aperture photometry, we also measure y a . The relative fluxes y a are 
depicted in Figure 7 (open star symbols). The Maidanak measurements in the first semester of 
2003 and the six Wise data of y a (open triangles; see here above) are tightly distributed around 
- 0.842 mag (solid line in Fig. 7). The rms fluctuation of the data is only of ~ 6 mmag (see the 
dashed lines in Fig. 7), which is consistent with the typical error of the measurements. The y a 
results in Fig. 7 suggest that both "a" and "b" are non-variable objects. Through 2003 (first and 
second semesters) we do not find any evidence in favour of variability of the "a-b" stars. 
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Fig. 9. Galaxy reconstruction obtained from selected Maidanak frames. Although the double 
quasar has been subtracted here, the positions of the components are labeled with two crosses: A 
is on the left and B is on the right, and the separation between both crosses is of 1'.' 1 - 1'.'2. The 
galaxy contours cover a region of 10" x 10". 

We show our global /?-band light curves of SBS 0909+532 in Figure 8. The open circles 
(Maidanak) and filled circles (Calar Alto) are the measurements of y A , whereas the open squares 
(Maidanak) and filled squares (Calar Alto) are the values of yb~ 0.45 mag. We have 3 1 points for 
the A component (circles) and 26 points for the B one (squares). The top panel of Fig. 8 contains 
the results in the winter-spring of 2003 and the bottom panel of Fig. 8 includes the results in 
the autumn of 2003. For each component we test the existence of a bias between the Calar Alto 
and Maidanak fluxes, e.g., f} A = y A (Calar Alto) -y A (Maidanak). Very small biases of /3 A = + 
15 mmag and /? B = - 30 mmag are found, and these corrections are taken into account to make 
the global records in Fig. 8. The biases are derived from the comparison between the Maidanak 
fluxes in a thirty day period (from day 2750 to day 2780) and the Calar Alto fluxes at equal or 
close dates (see the top panel of Fig. 8). 

To roughly estimate the contaminations from the direct PSF fitting technique, we take some 
of our best Maidanak images (in terms of seeing conditions, FWHM ~ 1 ") in the R band. A 
zoom-in of one of these best frames is shown in Fig. 2. Firstly, we combine the selected frames 
and derive a numerical model of the galaxy from a regularizing algorithm. To produce a more sta- 
ble reconstruction, the real galaxy profile is assumed to be close to the Sersic profile (Koptelova 
et al. 2005). Our deconvolution method differs only slightly from the former deconvolution tech- 
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niques by Magain, Courbin & Sohy (1998) and Burud et al. (1998). Figure 9 presents the galaxy 
reconstruction obtained from the stack of the /?-band selected frames. The box in Fig. 9 is 16'.'6 
on a side. The positions of the components are labeled with two crosses: A is on the left and B is 
on the right. The innermost contours are circular-elliptical rings, whereas the outermost contours 
show a less definite shape. Secondly, the selected frames are fitted to a photometric model that 
includes the galaxy brightness. Therefore, we are able to infer clean relative fluxes of A and B 
(without contamination by galaxy light) and to compare them with the contaminated ones (from 
direct PSF fitting). As result of the comparison, we report typical (averaged) contaminations of 
18.8 mmag and 4 mmag for the A and B components, respectively. These very weak contami- 
nations are in reasonable agreement with our preliminary considerations in the beginning of this 
section, and are taken into account in the measurement of the /?-band flux ratio in Section 5. 

4. Time delay 

To calculate the time delay between both components of SBS 0909+532, we use the R-bmd 
brightness records corresponding to the winter-spring of the year 2003. The /?-band records 
are more densely populated than the V-band ones. Moreover, the /?-band time coverage in the 
winter-spring of 2003 (about 120 days) is longer than the time coverage in the autumn of 2003 
(about 50 days). Thus we focus on the ,R-band data from day 2670 to day 2790, i.e., 22 points 
in the A component and 19 points in the B component (see the top panel of Figure 8). There 
are different number of points for component A and component B because we only consider 
fluxes with uncertainties below 40 mmag (see Section 3). As the B component is fainter, its 
photometric uncertainties are larger and the number of final data is smaller. The new light curves 
are characterized by a mean sampling rate of one point each six days. 

Once we have the data set, a suitable cross-correlation technique is required. Here we mainly 
use the x 2 minimization (e.g., Kundic et al. 1997) and the minimum dispersion (D 2 ) method 
(Pelt et al. 1994, 1996). However, although other techniques are probably less robust than the 
X 2 and D 2 ones (doing a first delay measurement, without a previous empirical determination), 
we also tentatively explore the modified cross-correlation function (MCCF) technique (Beskin 
& Oknyanskij 1995; Oknyanskij 1997). The MCCF combines properties of both standard cross- 
correlation functions: the CCF by Gaskell & Spark (1986) and the DCF by Edelson & Krolik 
(1988). We begin our analysis using the^ 2 method, which is based on a comparison between 
the light curve y A (or y B ) and the time shifted light curve y B (or y A ). For a given lag, one can 
find the magnitude offset that minimizes the x 1 difference. From a set of lags, it can be derived 
a set of minima (of^ 2 ), which permits to make a^ 2 spectrum: x 1 vs lag- The best solution of 
the delay is the lag corresponding to the minimum of the x 2 spectrum. In general, the shifted 
epochs t' B (or t' A ) do not coincide with the unchanged epochs t A (or t B ), so we estimate the values 
of yA(t' B ) (or y B (t' A )) by averaging the A (or B) fluxes within bins centered on times t' B (or t' A ) 
with a semiwidth a. To average in each bin, it is appropriate the use of weights depending on the 
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Fig. 10. x 1 spectra for a — 7 days (dashed line), a — 8 days (solid line) and a - 9 days (dotted 
line). The parameter a is the semi width of the bins in the A component, so the three spectra repre- 
sent different results of a cross-correlation with reasonable time-resolution (the mean sampling 
time is of 6 days/point). 
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Fig. 11. x 1 spectra for a — 9 days (dashed line), a = 10 days (solid line) and a = 11 days (dotted 
line), where a is the semiwidth of the bins in the B component. 

separation between the central time t' B (or t' A ) and the dates t& (or ts) in the bin. In principle, we 
concentrate in the interval [- 90, + 90] days, which includes the predicted negative delays (see 
Introduction) as well as a wide range of unlikely positive delays (positive delays are inconsistent 
with basic observations of the system). 

Firstly, the curve and the time shifted curve ys are compared with each other (using bins 
in the A component). In order to work with a reasonable time-resolution, we use a values less 
than or equal to two times the mean sampling time, i.e., a < 12 days. The x 2 value roughly 
grows with the size of the bin, and^- 2 ~ 1 for a = 7-9 days. For a = 7-9 days, there are best 
solutions Atba = + 46-48 days (x 2 = 0.97-0.98), and we show the corresponding spectra in 
Figure 10. We have drawn together the spectra for a = 7 days (dashed line), a — 8 days (solid 
line) and a — 9 days (dotted line). Apart from the main minima close to + 50 days, there are 
other secondary minima at negative and positive lags. In Fig. 10, two secondary minima seem to 
stay significant for all the bin sizes: the minima close to - 50 days and the probable edge effects 
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Fig. 12. Histograms from 1000 repetitions and the^f 2 minimization (bins in B and a = 10 days). 
Top panel: best solutions of the time delay. Bottom panel: best solutions of the magnitude offset 
(flux ratio). These distributions are consistent with a delay of about one and a half months and a 
time-delay-corrected flux ratio of about 0.59 mag. 

at + 80-90 days. We also compare the curve and the time shifted curve Va, using bins in the 
B component. For a = 10 days, we obtain a best solution Atba = - 44 days (x 2 = 1.15). Smaller 
and larger bins lead to solutions characterized by x 2 < 0.7 and^f 2 > 1.2, respectively. In Figure 
11, the solid line represents the spectrum for a = 10 days, while the dashed line represents the 
spectrum for a = 9 days and the dotted line traces the spectrum for a = 11 days. Main minima in 
the interval - 40-50 days appear in all these cases. Unfortunately, important signals at positive 
lags and probable border effects at + 80-90 days are again included in the complex spectra. The 
important structures at positive lags in Figs. 10-11 are probably caused by artifacts in the cross- 
correlation, so they have no physical origin, but are due to the 10/20-day gaps and the moderate 
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Fig. 13. Comparison between the shifted A light curve (circles) and the B light curve (squares). 
We use ax 2 minimization, bins in B and a — 10 days. Top panel: best solution (Atba = - 44 
days, AniBA = 0.59 mag). Bottom panel: solution for Atba = — 80 days (AmsA = 0.56 mag). 



variability of the components. Therefore, taking a — 10 days (bins in the B component) and a 
negative range [- 90, 0] days, we try to determine a pre-conditioned time delay. 

In order to derive uncertainties, we follow a simple approach. We make one repetition of the 
experiment by adding a random quantity to each original flux in the light curves. The random 
quantities are realizations of normal distributions around zero, with standard deviations equal 
to the errors of the fluxes. We can make a large number of repetitions, and thus, obtain a large 
number of At ba values. The true value will be included in the whole distribution of measured 
delays. From the^ 2 minimization (bins in B and a = 10 days) and 1000 repetitions, we obtain the 
histograms in Figure 12. Regarding the distributions in the top panel (delays) and bottom panel 
(flux ratios) of Fig. 12, the main features lead to measurements Atba = - 45 *J l days and Ahiba = 
0.590 + 0.014 mag (95% confidence intervals). We note that the main delay peak is asymmetric, 
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Fig. 14. Dispersion spectra for 5 = 1 days (dashed line), 5 = 9 days (solid line) and S = 11 days 
(dotted line). The parameter 6 is the decorrelation length associated with the D 2 A , spectra (Pelt et 
al. 1996). 

so 55% of the repetitions correspond to - 44-45 days, whereas 40% of the repetitions correspond 
to values < - 45 days. The secondary delay peak (around - 20 days) represents about 5% of 
the repetitions and is associated with the secondary minima in the negative region of Fig. 1 1 . 
Therefore, the distribution in the top panel of Fig. 12 permits a 95% estimation of the time delay 
of SBS 0909+532. 

In Figure 13 (top panel), the A light curve (circles) shifted by the optimal values of the 
time delay and the magnitude offset (time-delay-corrected flux ratio), and the unchanged B 
light curve (squares) are plotted. The cross-correlation using bins in the B component (a = 10 
days) indicates that the initial variations in the brightness of B reasonably agree with the final 
fluctuations in the brightness of A. The overlap for a delay of - 80 days (e.g., Saha et al. 2005) 
also appears in the bottom panel of Fig. 13. However, this last time delay is clearly rejected by 
the observations, since the;^ 2 value is larger than 10 (x 2 ~ 18). 

To confirm the results from the^ 2 minimization, we also use the dispersion spectra introduced 
by Pelt et al. (1994, 1996). The basic idea is a combination of yA and into one global record 
for every lag r and magnitude offset mo by taking all the values of yA as they are and shifting 
the values of yg - mo by r. For each t one can find the mo value that minimizes a dispersion 
estimate D 2 (t, mo), so a dispersion spectrum D 2 (t) can be made in a direct way. We focus on 
the D 2 A2 spectra that are called D 2 for simplicity (see Pelt et al. 1996 for details). This technique 
incorporates a decorrelation length (6), where 6 plays a role similar to that of a in the^ 2 method. 
Considering reasonable values of 6 (from 7 to 1 1 days, see here above), we are able to make some 
interesting spectra. In Figure 14 we have plotted together the spectra for 6 = 7 days (dashed line), 
6 = 9 days (solid line) and 6 = 1 1 days (dotted line). Although there are main minima in the 
interval - 40-50 days, there are also significant signals at positive lags and probable border 
effects at + 90 days. In the negative region of Fig. 14, a secondary minimum around - 70 days 
appears. Using 6 = 9 days and a negative range [- 90, 0] days, we carry out a second pre- 
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Fig. 15. Histograms from 1000 repetitions and the D 2 minimization (5 = 9 days). Top panel: 
best solutions of the time delay. Bottom panel: best solutions of the magnitude offset (flux ratio). 
These distributions are again consistent with an 1 .5-month delay and and a time-delay-corrected 
flux ratio of about 0.58-0.59 mag. 

conditioned measurement of the time delay. The uncertainties are deduced from 1000 repetitions 
of the experiment (see here above), and the relevant histograms are shown in Figure 15. While 
the top panel contains the distribution of delays, the bottom panel traces the distribution of flux 
ratios. Through the distributions in Fig. 15, we obtain that At ba = - 48 + J 6 days and Am BA = 0.585 
+ 0.020 mag (90% confidence interval). These D 2 results strengthen the conclusions from the^f 2 
technique. A marginal measurement (10% confidence interval) of At ba = - 67 ^ days and Am B A 
= 0.558 ^Q QQg mag is also possible. However, both this possibility and the^ 2 result of around - 
20 days are probably related to the presence of gaps and the absence of strong variability in the 
light curves. 
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A MCCF technique (Beskin & Oknyanskij 1995; Oknyanskij 1997) is also explored. The 
MCCF is a modification of the standard cross-correlation functions (CCF and DCF). When this 
MCFF is applied to our data in the lag interval [- 60, + 60] days, the maximum correlation 
coefficient (0.907) corresponds to a lag of - 45 days. This last result basically agrees with the^ 2 
and dispersion spectra in Figs. 1 1 and 14. 

5. Conclusions 

Nowadays several groups are trying to coordinate the rich but scattered research poten- 
tial in the field of gravitationally lensed quasar monitoring. The goals are to rational- 
ize the astronomical work and to catalyze big scientific collaborations so that the astro- 
physics community can get a significant progress in the understanding of the central engine 
in lensed quasars, the structure of the lensing galaxies and the physical properties of the 
Universe as a whole. Some examples about that are the Astrophysics Network for Galaxy 
LEnsing Studies (ANGLES, http://www.angles.eu.org/), the Cosmic Lens All-Sky Survey 
(CLASS, http://www.aoc.nrao.edu/~smyers/class.html) and the COSmological MOnitoring 
of GRAvItational Lenses (COSMOGRAIL, http://www.cosmograil.org/). The University of 
Cantabria group (Spain), three groups of the former Soviet Union (Institute of Astronomy of 
Kharkov National University, Ukraine, Sternberg Astronomical Institute, Russia, and Ulug Beg 
Astronomical Institute of Uzbek Academy of Science, Uzbekistan) and the Tel-Aviv University 
group (Israel) are also carrying out a series of initiatives to better exploit the recent individual 
monitoring campaigns as well as to solidify some future common project. In this paper we present 
the first collaborative programme on the variability of the double quasar SBS 0909+532A,B. The 
VR observations of the system and the field stars were made with three modern ground-based 
telescopes in the year 2003. 

The SBS 0909+532c star (N23210036195 in the GSC2.2 Catalogue) at (a, S) = (09: 12:53.59, 
+52:59:39.82) in J2000 coordinates is found to be variable, with two different levels of flux. 
The VR gap between the low-state and the high-state is of 80-100 mmag, and the low-state 
lasts about one month. In the high-state the star also seems to vary, but these small-amplitude 
variations are not so significant as the gap between states. We want to remark the variability of 
this nearby star ("c" star), and to encourage colleagues to follow-up its fluctuations and identify 
the kind of variable source. The "c" star cannot be used as the reference object (differential 
photometry), because it introduces a zero-lag global correlation between the light curves of the 
quasar components A and B. However, the "a-b" nearby stars are non-variable sources, and 
we choose the "b" star as the reference candle. On the other hand, the "si" and "s2" stars are 
relatively far objects, which were proposed as good references in a previous analysis (Nakos et 
al. 2003). However, the new /?-band light curve m S 2 - m s \ reveals the variability of one ("si" or 
"s2") or both stars. This variability could be either a very rare phenomenon or a consequence of 
doing more refined measurements (aperture or PSF fitting on several frames per night). We warn 
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about the possible problems with this pair of stars and think it merits more attention. The point- 
spread function (PSF) fitting methods permit to resolve the two components of the quasar and 
to derive the VR light curves of each component. These new VR light curves represent the first 
resolved brightness records of SBS 0909+532. Although the V-band curves are interesting, the 
/?-band records seem more reliable and are more densely populated. The /?-band curves show 
a moderate variability through 2003, and the observed fluctuations are promising for different 
kinds of future studies. 

To estimate the time delay between the components of SBS 0909+532, we use an 120- 
day piece of the /?-band brightness records, and x 2 and dispersion (D 2 ) techniques. The cross- 
correlation of the two light curves (A and B) leads to complex x 2 spectra. However, assuming 
that the quasar emission is observed first in B and afterwards in A, or in other words, Atba < 
(in agreement with basic observations of the system), 95% measurements Atba = - 45 *Jj days 
and AniBA = 0.590 + 0.014 mag are inferred from 1000 repetitions of the experiment (synthetic 
light curves based on the observed records). From the D 2 minimization (Pelt et al. 1996) and 
1000 repetitions, we also obtain 90% measurements Atba = - 48 + l days and ts.rn.BA = 0.585 + 

0. 020 mag. The D 2 uncertainties are derived under the already mentioned assumption that Atba 
is negative. There is a clear agreement between the results from both techniques, so a delay value 
of about one and a half months is strongly favoured. Our light curves rule out a delay close to 
three months, which has been claimed in a recent analysis (Saha et al. 2005). When we measure 
the time delay of the system, we simultaneously derive the time-delay-corrected flux ratio (at the 
same emission time) in the R band. This quantity, Aytlba = m^(f + Atba) - mA(f), is contaminated 
by light of the lens galaxy, and taking into account the weak contaminations of A and B (see the 
end of subsection 3.2), the totally corrected /?-band flux ratio is 0.575 + 0.014 mag. We remark 
that our final R flux ratio is in total agreement with the rough (uncorrected by the time delay and 
the contamination by galaxy light) measurement by Kochanek et al. (1997): 0.58 ± 0.01 mag. To 
properly determine a flux ratio, one must use clean fluxes at the same emission time, i.e., fluxes at 
different observation times and without contamination (Goicoechea, Gil-Merino & Ullan 2005). 
Only for particular cases (e.g., faint lens galaxy, short delay and moderate variability), it may be 
reasonable to use direct fluxes. 

In order to get a reasonably good value of x 2 , we do not need to introduce a time dependent 
magnitude offset or a complex iterative procedure (e.g., Burud et al. 2000; Hjorth et al. 2002), 

1. e., only a delay and a constant offset are fitted. This is a strong point of the analysis. The 
agreement between the results from different techniques is another strong point. However, the 
new measurements have some weak points that we want to comment here. The weakest point 
is the relatively poor overlap between the A and B records, when the A light curve is shifted 
by the best solutions of the time delay and the magnitude offset (e.g., see the top panel of Fig. 
13). Moreover, we carry out pre-conditioned measurements, since a negative interval [- 90, 0] 
days is considered in the estimation of uncertainties (component B leading component A). This 
second weak point is related to the presence of 10/20-day gaps and the moderate variability of 
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the components, which does not permit to fairly rule out positive delays. We nevertheless remark 
that the negative interval is in good agreement with the predictions by Lehar et al. (2000) and 
Saha et al. (2005), and we find^ 2 and D 2 minima around - 45 days when the observed data and 
both negative and positive lags are taken into account (see Figs. 1 1 and 14). Of course, as any 
another first determination of a time delay, the 1 .5-month value should be confirmed from future 
studies. 

Forty years ago, Refsdal (1964) suggested the possibility of determining the current expan- 
sion rate of the Universe (Hubble constant) and the masses of the galaxies from the time delays 
associated with extragalactic gravitational mirages. More recently, for a singular isothermal el- 
lipsoid (SIE), Koopmans, de Bruyn & Jackson (1998) found that the time delay can be cast in a 
very simple form, depending on basic cosmological parameters, redshifts and image positions. 
The relevant image positions are the positions with respect to the centre of the main lens galaxy, 
and the SIE delay is similar to the delay for a singular isothermal sphere (SIS). In principle, a 
singular density distribution is justified because a small core radius changes the time delay negli- 
gibly, and only a small core radius seems to be consistent with the absence of a faint central image 
(e.g., Kochanek 1996). Moreover, individual lenses and lens statistics are usually consistent with 
isothermal models (e.g., Witt, Mao & Keeton 2000 and references therein), so it is common to 
adopt an isothermal profile. Witt, Mao & Keeton (2000) showed that an external shear changes 
the simple SIS time delay in proportion to the shear strength. For two-image lenses that have a 
small shear and images at different distances from the centre of the lens, the shear should have 
a small effect on the time delay. Thus, when one has accurate measurements of image positions, 
redshifts and time delay, it is viable an accurate estimation of Ho (using complementary infor- 
mation on the matter/energy content of the Universe). 

Very recently, Kochanek (2002) also presented a new elegant approach to the subject. He 
modelled the surface density locally as a circular power law, with a mean surface density < k > 
in the annulus between the images. Expanding the time delay as a series in the ratio of the 
thickness of the annulus to its average radius, it is derived a delay that is proportional to the SIS 
time delay. The zero-order expansion term consists of the SIS delay and a multiplicative factor 
2(1- < k >). Kochanek also incorporated the quadrupoles of an internal shear (ellipsoid) and an 
external shear. However, for two-image lenses where the images lie on opposite sides of the lens, 
the delay depends little on the quadrupoles. This novel perspective is useful to infer < k > from 
observations of the lens system (time delay, image positions and redshifts) and complementary 
cosmological data (expansion and matter/energy content of the Universe). 

For SBS 0909+532, although the redshifts are very accurately known and the time delay is 
now tightly constrained (or at least there is a first accurate estimation to be independently con- 
firmed), the inaccurate position of the main lens galaxy does not permit an accurate measurement 
the cosmic expansion rate and the surface density of the main deflector. We have Hq oc q\ - 0^ 
and 1- < k >oc (0^ - fl^) -1 , where 6 A and 9 B are the image angular positions with respect to the 
centre of the main lens galaxy. On the other hand, using the astrometry in Table 3 of Lehar et al. 
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(2000), it is easy to obtain & 2 B - Q 2 A - 0.4 + 0.2. Thus we conclude that the accuracy in 0^ - 6 2 A is 
only 50%, indicating the necessity of new accurate astrometry of SBS 0909+532. 
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